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The process of jet-gap-jet (JGJ) production is briefly described. The 
JGJ scattering amplitude parametrisation is discussed. On the basis of 
full amplitude calculations, the parametrisation formulas for the leading 
logarithmic (LL) and next-to-leading logarithmic (NLL) approximations 
are obtained. For each case a sum over all conformal spins is considered. 
The obtained agreement is better than 0.25% for LL and 1% for NLL. 
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1. Introduction 

Hard diffractive processes have been an important part of the studies 
performed in high energy physics since their discovery in the UA8 experi¬ 
ment [T]. The data collected by HERA and Tevatron detectors allowed to 
deepen these studies. Nevertheless, now, in the LHC era, many questions 
still remain open. 

The definition of diffraction is usually connected with an exchange of a 
colourless object. In the case of the electromagnetic fields the exchange is 
mediated by a photon, whereas in the case of strong interactions it is related 
to the Pomeron. A colourless exchange leads to one of the most prominent 
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features of diffraction - the presence of a large rapidity gap. In theory, a gap 
is a space interval in rapidity in which no particles are produced. Moreover, 
since the colourless exchange does not influence the quantum numbers, the 
state of interacting objects is preserved. 

In particular, jet events could be created as a result of the colourless 
exchange between interacting gluons (see Fig. [T|). Such a process will be 
hereafter called the jet-gap-jet (JGJ) production. 



Pomeron 



Fig. 1. Diagram of the jet-gap-jet production. 


The presence of a hard scale in the jet-gap-jet events makes it possi¬ 
ble to understand these events with perturbative methods. However, after 
many theoretical investigations, there is still no consensus on what is the 
relevant QCD mechanism. When the rapidity gap is sufficiently large, the 
perturbative QCD description of the jet-gap-jet events is usually performed 
in terms of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) Pomeron [2]. As 
was shown in Ref. [3], this model can be tested experimentally by studying 
the behaviour of the ratio of dijet events with a gap to all dijet events as a 
function of the leading jet transverse momentum or the interval in rapidity 
between the two leading jets. 

It is worth mentioning that the jet-gap-jet events can also be produced in 
single diffractive and double Pomeron exchange topologies [3]. Such events, 
so far never observed experimentally, are expected to be measured at the 

LHC. 

In order to account for the experimental effects, the detector simulation 
is usually performed. The events used in such simulations should have prop¬ 
erties similar to the ones that are to be measured. This requires a use of the 
Monte Carlo (MC) generators. The process of the jet-gap-jet production was 
implemented in several such tools, e.g. Herwig [5] or FPMC [6]. In order to 
speed up the calculations, the cross-section formulas were parametrised [3]. 
These equations were fitted assuming that the transverse momentum of the 
leading jet did not exceed 120 CeV. However, as this assumption is no longer 
valid at the energies available at the LHC, a new parametrisation has to be 
carried out. In this paper we present new fit formulas applicable for the 
LHC energies. 
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2. Jet-Gap-Jet Formalism 

The parton-level cross-section for the jet-gap-jet production can be cal¬ 
culated as [3]: 


dSfjPP^XJJY 

dxidx2dp^ 


■ feff{xi,PT) ■ feff{x2,PT) 


da99^99 

dp'^ 


where 5^ is the gap survival probability, feffixi, 2 ,PT) are the effective Par- 
ton Distribution Functions (PDFs): 


feff{xi,2,PT) = 9 {x,Pt) + {q{x,PT) + q{x,PT)) , 


where Nc = 3 is the number of colours, Cp - the QCD colour factor and 
g, q and q are the gluon, quark, antiquark distribution functions in the 
interacting hadrons. The cross-section is given by: 


da99^99 
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where Xeff ~ the effective BFKL kernel [3], a‘^{p‘^) = TTa1{p‘^)/Nc is the 
running coupling constant and p denotes the conformal spin. The complex 
integral runs along an imaginary axis from ^ —ioo to ^ -|-ioo. The necessity 
of summing up all the conformal spins was demonstrated in [7] . The BFKL 
kernel was calculated in the leading logarithmic (LL) [2] and next-to-leading 
logarithmic (NLL) [8j approximations. 


3. Jet-Gap-Jet Amplitude Parametrisation 

In principle, the JGJ amplitude can be directly implemented into the 
Monte Carlo generator. Unfortunately, due to the complexity of Eq. [H the 
computation time is quite long. Since the usual number of events for a 
physics analysis is of the order of 10®, the parametrisation procedure was 
postulated to speed up the generation process [3]. This parametrisation 
was done for the phase-space expected to be measurable using the Tevatron 
data [9]: 

• transverse momentum of jets: 20 < pp < 120 GeV, 
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• pseudorapidity distance between jets: 0 < Ar/ < 10. 

The conformal spins were summed from p = —10 to p = 10, which is 
sufficient for the considered jet px range [3]. The parametrisation procedure 
resulted in the speed-up factors of 10 and 1000 for LL and NLL, respectively. 
In the following, we redo the parametrisation for the values expected to be 
observed at the LHC |10j : the range of jet transverse momentum extended 
to 1 TeV and the sum over conformal spins from -50 to 50. 

3.1. Leading Logarithmic Approximation 

The leading logarithmic approximation is known to be insufficient as the 
next-to-leading BFKL corrections are expected to be large [S]. However, for 
the completeness, we discuss it below. Denoting z{p‘^) = a^{p‘^) ■ Arj/l, the 
LL cross-section can be parametrised as: 

• = A • [^ -|- exp(H + C- z + D-z^ + E-z^ + F- z^)], for p = 
0, 

• A'^\^^{z) = N ■ [A + B ■ z + exp(C' + D- z + E-z‘^ + F- , for sum 
over conformal spins, 

_ 2 

where the normalisation constant is equal to A = ^ with fixed to 0.17. 

The shape of the full amplitude as a function of pseudorapidity difference 
is shown in Fig. [2] (top). In the bottom, the comparison between the full 
amplitude calculations and the parametrisation results (fit) is presented. 
The obtained fit parameters are listed in Table [H For both considered cases 
the differences are well below 2.5%o for the whole pseudorapidity range. 

Table 1. Fit parameters for leading logarithmic amplitude parametrisation for 
p = 0 and sum over conformal spins (all p). 


Parameter 

p = 0 

all p 

A 

0.452 ± 0.023 

-2.032 ± 0.022 

B 

2.2162 ± 0.0022 

-1.135 ± 0.057 

C 

-6.436 ± 0.029 

6.18035 ± 0.00012 

D 

21.16 ± 0.12 

-14.3093 ± 0.0032 

E 

-16.46 ± 0.14 

20.650 ± 0.011 

F 

5.586 ± 0.056 

-6.4983 ± 0.0096 
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Fig. 2. Top: the shape of the full amplitude as a function of pseudorapidity differ¬ 
ence. Bottom: comparison of the leading logarithmic full amplitude calculations 
and parametrisation results (fit). Black, solid line is for the sum over all conformal 
spins and the red, dashed one is for p = 0. 


3.2. Next-to-Leading Logarithmic Approximation 

The parametrisation of NLL amplitude is more complex as it depends on 
both: jets transverse momenta (identical for both jets in collinear approach) 
and their pseudarapidity distance. In order to obtain the parametrisation 
formulas, the one-dimensional problem of finding the general dependence on 
the jet transverse momentum was addressed first. This formula was found 
to be: 


^AfLL(PT, Ary = fixed) 


N ■ 


A{Ar]) ■ -h C{Ar]) ■ 


( 2 ) 


The exemplary results obtained for three different pseudorapidity distances 
(Ary = 0.5, Ary = 5.5 and Ary = 9.5) are shown in Fig. [3l In the top 
part of this figure the shape of the distribution is presented, whereas in the 
bottom the ratios of full amplitude calculation to fit are shown. The overall 
agreement is well below 1%, which means that Eq. [2] works for the whole 
considered Ary range. 

The parametrisation formulas for both cases (p = 0 and the sum over 
all conformal spins) were found to have the following forms: 

• A{z) = ao + ai ■ z + exp(a 2 + ■ z + ■ z'^ + ■ z^), 

• B{z) = bo + bi- z, 

• C{z) = Co + Ci - Z + exp(c 2 + C 3 ■ Z + C 4 ■ z"^ + C 5 ■ z^), 
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Fig. 3. Top: the shape of the full amplitude as a function of jet transverse momen¬ 
tum. Bottom: comparison between next-to-leading logarithmic full amplitude 
calculations and parametrisation results (fit) for p = 0 (left) and sum over all con¬ 
formal spins (right). Black rectangles are for At] = 0.5, red circles for Ar] = 5.5 
and blue triangles for Ary = 9.5. 


• D{z) = do + di ■ z + d 2 ■ z'^ + ds ■ . 

The shape of the full amplitude as a function of pseudorapidity difference 
for the NLL is shown in Fig. 0] (top). In the bottom, the comparison of the 
full amplitude calculations and the parametrisation results is shown. The 
obtained fit parameters are listed in Table [21 For both considered cases the 
differences are below 1% for the whole considered pseudorapidity range. 


4. Summary 

The process of jet-gap-jet production is very interesting to be studied. 
The investigation of its properties would allow us not only to measure the 
cross-section, but also to test the BFKL model. 

The previous parametrisation, constructed to describe the Tevatron data, 
was extended to include phase-space regions available at the LHC energies. 
The jet transverse momentum was assumed to be within the 20 GeV to 
1 TeV range and the conformal spins were summed from -50 to 50. The ob¬ 
tained agreement between the full amplitude calculations and the parametri¬ 
sation results was found to be better than 0.25% for leading logarithmic and 
1% for next-to-leading logarithmic approximations. 
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Fig. 4. Top: the shape of the full amplitude as a function of jet transverse momen¬ 
tum. Bottom: Comparison between next-to-leading logarithmic full amplitude 
calculations and parametrisation results (fit) for p = 0 (left) and sum over all con¬ 
formal spins (right) as a function of presudorapidity difference. Solid black lines are 
for jet transverse momentum of pr = 20 GeV, red dashed ones for pr = 320 GeV, 
green fine-dashed for px = 620 GeV and blue dotted for px = 920 GeV. 
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